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Abstract 

Based on SU(3) flavor symmetry, many of the quark-graph amplitudes in two-body nonleptonic 
decays of charmed mesons can be extracted from experiment, which enable us to see the relevance 
and importance of weak annihilation topologies and to determine the complex parameters a\ and 
a2 to test the factorization approach. It is found that a2 in D — > K*n and D — > Kp can be different 
by a factor of 2, indicating that nonfactorizable corrections to the latter are far more important 
than the former. The relative phase between a\ and 02 is about 150°. Weak annihilation topologies 
induced by nearby resonances via final-state rescattering can be described in a model-independent 
manner. Although the W-exchange contribution in D —* PP decays is dominated by resonant 
final-state interactions (FSIs), its amplitude in VP decays (V: vector meson, P: pseudoscalar 
meson) receives little contributions from FSIs in the quark- antiquark resonance formation. As a 
consequence, the sign flip of the PF-exchange amplitude in D — > K*tt and Kp decays, which is 
needed to explain the relatively real decay amplitudes of D — > Kp, remains unexplained. SU(3) 
symmetry is badly broken in some Cabibbo-suppressed modes and it can be accounted for by the 
accumulation of some modest SU(3) violation in individual quark-graph amplitudes. 
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I. INTRODUCTION 



The hadronic decays of charmed mesons and related physics have been studied extensively 
in past 25 years and a lot of progress has been made. The charm lifetimes, e.g., t(D^) and 
t(S+), the D° — D° mixing and the Dalitz plot analyses of three-body charm decays are 
some of the main topics that are currently being studied. Many new results are expected 
soon from the dedicated experiments conducted at CLEO, E791, FOCUS, SELEX and the 
B factories BaBar and Belle. 

Contrary to the experimental progress, the theoretical advancement is relatively slow. It 
is known that the conventional naive factorization approach fails to describe color-suppressed 
(class-II) decay modes. Empirically, it was learned in the 1980s that if the Fierz-transformed 
terms characterized by 1/N C are dropped, the discrepancy between theory and experiment is 
greatly improved JT[. This leads to the so-called large- N c approach for describing hadronic 
D decays 0. Theoretically, explicit calculations based on the QCD sum-rule analysis @ 
indicate that the Fierz terms are indeed largely compensated by the nonfactorizable correc- 
tions. Due to the success of the 1/N C approach to charmed meson decays, it was widely 
believed in the eighties that it applies equally well to the weak hadronic decays of bottom 
mesons. However, a generalization of the large-iV c approach or the sum rule analysis |4j] to 
hadronic B decays leads to some predictions in contradiction to experiment, namely, the 
destructive interference in the class-Ill modes B~ — > D 0( -*^ir~ is not borne out by the data. 
In the heavy quark limit, nonfactorizable corrections to nonleptonic decays are calculable 
due to the suppression of power corrections. Unfortunately, the charmed quark is not heavy 
enough to apply the QCD factorization approach H or pQCD in a reliable manner. 

Moreover, the importance of weak annihilation contributions, namely, VF-exchange and 
^-annihilation, is controversial. In practical calculations, it is customary to argue that 
they are negligible based on the helicity suppression argument. Although the observation 
of D° — > ~K°<j) in the middle 1980s seems to give the first clean evidence of lU-exchange, 
it was claimed in || that rescattering effects required by unitarity can produce the same 
reaction even in the absence of the VT-exchange process. Then it was shown in that this 
rescattering diagram belongs to the generic iy-exchange topology. It has been stressed in 
|7]] that even in D — > Kir decays, the ^-exchange contribution is sizable. 

It has been established sometime ago that a least model-independent analysis of heavy 
meson decays can be carried out in the so-called quark-diagram approach. In this diagram- 
matic scenario, all two-body nonleptonic weak decays of heavy mesons can be expressed in 
terms of six distinct quark diagrams T, the color-allowed external W-emission tree 

diagram; C, the color-suppressed internal H^-emission diagram; S, the PU-exchange diagram; 
A, the H^-annihilation diagram; V, the horizontal H^-loop diagram ; and V, the vertical 



* Historically, the quark-graph amplitudes T, C, £, A are originally denoted by A, £>, C, V, respec- 
tively lil. 
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M^-loop diagram. (The one-gluon exchange approximation of the V graph is the so-called 
"penguin diagram".) It should be stressed that these quark diagrams are classified according 
to the topologies of weak interactions with all strong interaction effects included and hence 
they are not Feynman graphs. All quark graphs used in this approach are topological and 
meant to have all the strong interactions included, i.e. gluon lines are included in all possi- 
ble ways. Therefore, topological graphs can provide information on final-state interactions 
(FSIs). 

Based on SU(3) flavor symmetry, this model- independent analysis enables us to extract 
the topological quark-graph amplitudes and see the relative importance of different under- 
lying decay mechanisms. The quark-diagram scheme, in addition to be helpful in organizing 
the theoretical calculations, can be used to analyze the experimental data directly. When 
enough measurements are made with sufficient accuracy, we can find out the values of each 
quark-diagram amplitude from experiment and compare to theoretical results, especially 
checking whether there are any final-state interactions or whether the weak annihilations 
can be ignored as often asserted in the literature. 

Recently, Rosner |nj has determined the diagrammatic amplitudes from the measured 
Cabibbo-allowed two-body D decays. There are several important observations one can learn 
from this analysis. First, the weak annihilation (W-exchange or W- annihilation) amplitude 
is sizable with a large phase relative to the tree amplitude. Second, the three D — > Kp 
amplitudes are observed to be relatively real, in sharp contrast to Kir and K*ir cases. It was 
argued in [|10j that the ^-exchange amplitude has to flip its sign from K*ir to Kp in order to 
explain why the three decay amplitudes of D — > Kp are in phase with one another. Third, 
the color-suppressed amplitude C has a nontrivial phase relative to the tree amplitude T . 
As we shall see, the appearance of nontrivial relative phases between various quark-graph 
amplitudes implies the relevance and importance of inelastic final-state interactions. 

The purpose of this work is twofold: First, we will utilize the reduced quark-graph 
amplitudes extracted from the data to determine the complex parameters a\ and a<i appearing 
in the factorization approach. This enables us to test the factorization hypothesis and see 
how important the nonfactorizable correction is. Second, we will study weak annihilations 
induced from nearby quark-antiquark intermediate states. This allows us to explore the 
effect of inelastic FSIs and see if the sign of the ^-exchange topology in Cabibbo-allowed 
D — > VP decays is governed by nearby resonances. 

The layout of the present paper is as follows. In Sec. II we first discuss the quark-diagram 
amplitudes and then extract the parameters a% and 02. The diagrammatic amplitudes for 
Cabibbo-suppressed decay modes and SU(3) violation are addressed. The weak annihilation 
induced from final-state rescattering in resonance formation is studied in Sec. III. Its impli- 
cation and importance for explaining some D decay modes are discussed. Sec. IV is devoted 
to exploring the color-suppressed amplitude and its phase. We then compare the present 
study with B decays in Sec. V and give conclusions in Sec. VI. 
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II. DIAGRAMMATIC APPROACH 



A. Quark-graph amplitudes 

Based on SU(3) flavor symmetry, the quark-graph amplitudes for the Cabibbo-allowed 
decays of charmed mesons are listed in Table I. (For a complete list of the quark-graph am- 
plitudes for Cabibbo singly and doubly suppressed modes, see f7j.)[] Note that the selection 
rule for a vanishing Df — > 7r + 7r° follows from the isospin transformation properties of the 
weak Hamiltonian and isospin invariance of strong interactions and hence it is unaffected by 



SU(3) breaking or final-state interactions |14[]. For final states involving 77 or 7/ it is more 



convenient to consider the flavor mixing of rj q and r] s defined by 

rj q = -^( uU + dd), r] s = ss, (2.1) 

in analog to the wave functions of to and <fi in ideal mixing. The wave functions of the 77 and 
77' are given by 



rj \ I cos <p — sin (p \ I 7] q 
rf J \ sin (f) cos (j) J \Vs 



(2.2) 



where <fi = 9 + arctanv^2, and 9 is the 77 — 77' mixing angle in the octet-singlet basis 

/ r] \ ( cos 9 — sin 9 \ / r] 8 \ 
\ rj ' J \ sin 9 cos 9 J \Vo J 

The D — > Mr) and M77' amplitudes have the expressions 



(2.3) 



A{D -> M77) = A(D -> Mt7 8 ) cos0 - A(D -> Mr/ ) sin0, 
A(D -> Mr]') = A(D -> Mr] 8 ) sin + A(D -> M/7 ) cos 0. (2.4) 

Note that the D — > i^7r amplitudes satisfy the isospin triangle relation 

A(D + -> ^°7T + ) = A(Z)° -> ^-7T + ) + ^(D -> K°7T ) (2.5) 

and likewise for D — > i^*7r and £) — > Kp amplitudes. Now since all three sides of the isospin 
triangle are measured, we are able to determine the relative phases between various decay 
amplitudes. From the measured decay rates Hl5|, we find (only the central values are quoted) 



^For charm decays involving an SU(3) singlet in the final product, e.g., D° — > K°<p, ~K°lo, K°t]o, 
there exist additional hairpin diagrams in which a quark-antiquark pair is created from vacuum 
to form a color- and fiavor-singlet final-state meson |lljl2|. There are four different types of 
disconnected hairpin diagrams: Sh, Ah,Vh, T^h corresponding to the quark graphs £, A, V, V (for 
details, see [ll]]). Here we will omit the contributions from the hairpin diagrams, though they seem 
to play some role in Df —> p + r)' [O] and B — > Krj decays |Q. 



4 



TABLE I. Quark-graph amplitudes for Cabibbo-allowed decays of charmed mesons. For re- 
duced amplitudes T and C in D — ► VP decays, the subscript P (V) implies a pseudoscalar (vector) 
meson which contains the spectator quark of the charmed meson. For £ and A amplitudes with 
the final state (71(72, the subscript P (V) denotes a pseudoscalar (vector) meson which contains the 
antiquark 



D + -► K°ir+ T + C 


D+ -► K* 7T+ T V + C P 




_3rV j= 2 (c-£) 

- f 


D°^K*-tt+ T V + £ P 

-> K r) s £y 


D°^K-p+ T P + £ V 

^K°p° j= 2 (Cv-£v) 
—> K°lo j=(C v + £ v ) 
-► £ p 


£>+ -► K u if + C + A 
-► vr+vr 


Df -► iT* U K+ C P + ^4p 

^P+vr ^(--4p + A0 

-> p + ?? s Tp 


£>+ -► if U K*+ Cy + -4y 

-^V j=(A P -Av) 

-^TT+LO -t=(Ap + A V ) 
— > 7T + </> 7\/ 



^P>0- > i?°7r0,D0^^-7r+ ~~ 30 °' ^D+^X° 7 r+,D0^^-7r+ ~~ 80 °' 

where we have used the relation, for example, 



B(D° -> K-7T+) + 2B(D° -> TfV) - I & ) B(D + -> tfV) 

cos 5 r o = _ n (2.7) 

1,1 2^/B(D^ K-n+)^2B(D° -> ZV) 

to extract the phases. Therefore, the isospin triangle formed by amplitudes is dramat- 
ically different from the one constructed by i\V or K tt. This triangle is almost flat with 
zero area, indicating that the three Kp amplitudes are relatively real. This also can be seen 
from the isospin analysis: 

A(D° - K-n + ) = s^-A 1/2 + ^3/2, 

A(D° -> tfV) = -^|a 1/2 + (2.8) 
^ TfV) = VsA 3/2 , 

with = |Aj|e l<5i . It turns out that the isospin phase difference is consistent with zero for 
D — > ~Kp (see Table II). 
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TABLE II. Isospin amplitudes and phase differences for D — ► Kit, K*ir, Kp decays. 





D ^Ktt 




D^Kp 


1^1/2/^3/2 


3.83 ±0.27 


5.61 ±0.35 


3.59 ± 0.69 


1^1/2 ~ ^3/2 


(90 ±6)° 


(104 ± 13)° 


< 27° 



The reduced quark-graph amplitudes T, C, 8 have been extracted from Cabibbo-allowed 
decays by Rosner |L0| with the results:f] 

T = 2.69 x 10~ 6 GeV, C = (1.96 ± 0.14) e^ 152 " x 10~ 6 GeV, 

8 = (1.60 ± 0.13) e ill4 ° x 10~ 6 GeV, A = 1.10 e~ i70 ° x 10" 6 GeV (2.9) 

from D — ► i^7T, ^°?7, -ft' //' decays, and 

T y = (1.78 ± 0.22) x 10~ 6 , C P = 1.48 e" il52 ° x TO -6 , 8 P = (1.39 ± 0.08) e m ° x 10" 6 (2.10) 



from D — > if 7r, K <p and — ► n + (f). Without loss of generality, T has been chosen to be 
real. The amplitudes T, C,8,A have dimensions of energy as they are related to the decay 
rate via 

V{D^PP) = -^ r \A\\ (2.11) 

with p c being the cm. momentum. In contrast, the reduced amplitudes T P y,Cpy,8py are 
dimensionless as they are extracted from the relation 

3 

T(D -> VP) = -^—\A\ 2 . (2.12) 

8-Krriy 

Note that our convention for D — > V^P amplitudes is the same as that in J7J but different 
from those of Rosner [[10 1 ; this allows one to compare the theoretical calculations of VP 
amplitudes directly with the quark-graph amplitudes extracted from experiment. It should 
be stressed that the solutions given above are not unique. For example, a small T amplitude 
with the magnitude of 1.1 x 10 -6 GeV is also allowed. However, it is not favored by the 
factorization approach |lu]|7]|. The original analysis by Rosner is based on the 77 — rf mixing 
angle 9 = —19.5° (or = 35.2°) from which the 77 and rf wave functions have simple 
expressions |L6 



J The value of the ^-annihilation amplitude A is slightly different from that given in 
have used the D+ lifetime: t(D+) = (0.496±g;gjg) x 10~ 12 s 
difference \5ae\ quoted in Table IV of |L0| is too large by 10°. 



] as we 

It seems to us that the phase 
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V = ~^(^9 ~ Vs) = -^{uu + dd- ss), 

rj' = ~j={v q + V2r] s ) = -^{uu + dd + 2ss). (2.13) 
a/3 v6 



From Table I, it is easily seen that A(D° — > K rj) — C/y/3. Hence, the diagrammatic 
amplitude C is ready to be determined once this mode is measured. A phenomenological 
analysis of many different experimental processes indicates 9 = —15.4° or <fi = 39.3° [jl7|| . 
However, we find that the above diagrammatical amplitudes Q2.9|) and ( |2.10|) describe the 
observed rates well. 

In order to extract the quark-graph amplitudes Tp, Cy and Ey from D — > Kp, ~K*r], 
K*rf , K°u decays, some assumptions have to be made. As noted in passing, the most 
prominent feature of the D — > Kp data is that all their three decay amplitudes are almost in 
phase with one another. Therefore, the quark-graph amplitudes have to be aligned in such 
a way to render the resultant various decay amplitudes of D — > Kp parallel or antiparallel. 
There exist three possibilities: (i) Tp,Cy and £y are relatively real, (ii) the amplitudes Tp 
and Cy possess a relative phase of order —150° as in Kit and K*ir cases, but £y ~ —£p, and 
(iii) Ey is close to Ep, but the relative phase between Cy and Tp changes a sign. The first 
possibility was first pointed out by Close and Lipkin fll8]. It turns out that the four data of 



D — > Kp and D° — > K°u can be fit by setting the three quark-graph amplitudes to be 

(i) T P = 1.20 x 10" 6 , C v = 0.21 x 10" 6 , Ey = 1.56 x KT 6 . (2.14) 

As stressed by Close and Lipkin, this fit implies that Ey > Tp ^> Cy and that the interference 
in the decay D + — > K p + is constructive, contrary to the case of D + — > K n + and K n + . 
However, this fit will be discarded for the following reason. Since Tp < Ty and Cy Cp, the 
branching ratios for — > K K* + and D+ — > p + (r),r)') will become quite small. Neglecting 
iy-annihilation amplitudes Ay and Ap for the moment, the fit ( [2.14 ) leads to the predictions 



B(D+ -> TCK* + ) = 5.3 x 10- 4 , B(D+ p+rj) = 1.1%, B(D+ p+rj') = 0.45%, (2.15) 

which are too small compared to the corresponding experimental results |15|| : (4.3 ± 
1.4)%, (10.8 ± 3.1)%, (10.1 ± 2.8)%. As shown below, A P and Ay are constrained by 
the measurements of — > 7r + p° and tt + uj and they are small in magnitude. Within the 
allowed regions of Ay and Ap constrained by Eqs. ( |2.17| ) and ( |2.18| ), the predicted branching 
ratios for aforementioned three modes are still too small, especially for the K K* + decay. 
For example, it is found that B(D+ -> K°K*+) = 0.35% for Ay = 3.3 x lO" 7 . Fit (i) is 
also unnatural in the sense that FSIs will generally induce relative phases between various 
quark-graph amplitudes. 

The second scenario is considered by Rosner [TTJ] based on the argument that if the W- 
exchange amplitude is dominated by quark- ant iquark intermediate states, then a sign flip 
of Ey relative to Ep will be a consequence of charge-conjugation invariance. As stressed by 
Rosner, the presence of large final-state phases in the Kp case is masked by the cancellation 
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between Cy and Sy. This accidental cancellation arises in Kp decays but not in Kir and 
K*ir decays. In the following we list other two fits for the quark-graph amplitudes of Tp, Cy 
and Cp\ 

(ii) T P = 2.24 x 10~ 6 , C v = (1.07 ± 0.18) e - a48 ° x 1(T 6 , Sy = 1.20 e" l72 ° x 10~ 6 , 
(iii) T P = 2.24 x 10~ 6 , C v = (1.07 ±0.18)e il48 ° x 10" 6 , £ v = 1.20 e i72 ° x HT 6 , (2.16) 



where fit (ii) was first obtained by |]10|| . It is easily seen that although both fits give the same 
decay rates for D — > Kp, they yield different predictions for D° — > K* rj. For = 39.3°, 
we find B(D —> K rf) = 1.5% for fit (ii) and 0.78% for fit (iii), while the experimental 
branching ratio is (1.9±0.5)% fE5|. Hence, it appears that fit (ii) for quark-graph amplitudes 



is preferred. Moreover, as we shall see in Sec. IV, a model calculation of inelastic final- 
state rescattering indicates that the imaginary part of the color-suppressed amplitude Cy 
is negative, in accord with fit (ii). However, it will be shown later (Sec. III.C) that the 
W-exchange amplitudes Sp and £y are not dominated by resonant FSIs and hence the sign 
flip of Sy from Sp remains unexplained. Therefore, fit (iii) for quark-graph amplitudes Cy 
and Sy is not entirely ruled out. It is worth mentioning that, contrary to fit (i), fit (ii) or (iii) 
has relatively large tree and color-suppressed amplitudes. As a consequence, the interference 
occurring in D + — > K p + has to be destructive in order to accommodate the data. 

For the H^-annihilation amplitudes Ap and Ay, the measurement of — > tt + uj |I5| 
leads to 

\A P + Ay\ = (4.5 ± 1.0) x 10~ 7 , (2.17) 



while the amplitude \Ap — Ay\ can be extracted from the recent E791 experiment |T9| 
T(D+ -> P °tt + )/T{D+ -> n+n+n-) = (5.8 ± 2.3 ± 3.7)%, though it does not have enough 
statistic significance. The result is 

\A P -Av\ = (2.0 ±1.2) x 10~ 7 , (2.18) 



where use of B(D+ — ► -k + it + -k~) = (1.0 ± 0.4)% [jl5[ has been made. It will be shown in 
Sec. III.B that Ap — Ay receives dominant contributions from resonance- induced FSIs. Eqs. 
( p. 17] ) and ( |2.18|) suggest that the phase difference between Ap and Ay is less than 90° and 



that the magnitude of Ap or Ay is smaller than Sp and Sy, contrary to the PP case. The 
suppression of iy-annihilation will be explained in Sec. III. 
Without VT-annihilation, fit (ii) leads to 

B(Df -> K°K* + ) = 1.4%, B(D+ p + ri) = 3.9%, B(D+ -»• p+rf) = 1.6%. (2.19) 

In the presence of IV-annihilation contributions, the decay amplitudes of Df — ► p + r]^ read 
[see Table I and Eq. (gj)] 

A(Dj -> p + r]) = T P sin <p \={A P + Ay) cos0, 

v2 

A(D+ -> pV) = Tp cos (f>+ -^=(A P + Ay) sin 4>. (2.20) 

v2 
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It is obvious that p + r] and p + r]' cannot be accommodated simultaneously because if W- 
annihilation contributes constructively to the former, it will contribute destructively to the 



latter, and vice versa. This issue has been discussed in |20|-f23l and it is generally believed 
that the large discrepancy between theory and experiment means that there is an additional 
contribution to p + r]' owing to the special character of the rj . As an illustration, we find 

B(D+ -> K°K* + ) = 2.2%, B(Dj -> p + v ) = 5.4%, B(D+ -»• p+rf) = 1.2%, (2.21) 

for Ap + Ay = —4.5 x 1CT 7 and Ay = —3.3 x 1CT 7 . The smallness of p + rj' from tree 
and ^-annihilation contributions may indicate the relevance and importance of the hairpin 
diagrams for the p + i]' decay. For example, an enhancement mechanism has been suggested in 
pl| that a cs pair annihilates into a W + and two gluons, then the two gluons will hadronize 
mostly into rf . The other possibility is that the gluonic component of the rf , which can be 
identified with the physical state e.g. the gluonium, couples to two gluons directly. 

There are several important observations of the above extracted reduced quark-graph 
amplitudes: (i) The W -exchange or ^-annihilation contribution is in general comparable 
to the internal P^-emission and hence cannot be neglected, as stressed in 0. The weak 
annihilation amplitude has a phase of order 90° relative to T and this is suggestive of 
the importance of resonant contributions to weak annihilations in D decays, (ii) The W- 
annihilation A and jy-exchange amplitudes £ have opposite signs, (iii) The color-suppressed 
internal P^-emission amplitude C has a phase ~ 150° relative to T for all Cabibbo-allowed 
D decays. In the factorization approach, the relative phase is 180°. (iv) The W-exchange 
amplitude in K*ir and in Kp has an opposite sign. As stressed by Rosner, this sign difference 
is very crucial to explain why D — > Kp amplitudes are relatively real, but not the case for 
D — > K*tt and D — > Kir decays. 

It was conjectured in []10| that (i) the opposite sign between Sp and Sy arises from the 
fact that they are dominated by the quark- ant iquark intermediate states which have equal 
and opposite couplings to K*~tt + and K~p + by charge-conjugation invariance, and (ii) the 
relative phase between C and T comes from inelastic final-state rescattering. We will come 
to these issues in Sees. Ill and IV. 



B. Parameters a\ and 0,2 

In terms of the factorized hadronic matrix elements, one can define a% and 02 by 
T = U< - m 2 K )F DK (ml), 

C = ^=V ud V; s a 2 (K7r) f K (m 2 D - m^^K), 
T v = ^V ud V;M~K*7c)2f 7T m K *A° K \ml), 

C P = ^V ud V: s a 2 (K\) 2f K *m K *F^{m 2 K *), (2.22) 
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T P = ^V ud V;MK P )2f p m p F 1 DK (ml 



^V ud V;MKp) 2f K mX P «), 



where we have followed j24j for the definition of form factors. Factorization implies a uni- 



versal, process-independent a\ and a 2 , for example, a 2 {Kp) = a 2 (K it) = a 2 (Kir). 

In order to extract the values of a\ and a 2 we consider two distinct form factors models: 
the Bauer-Stech-Wirbel (BSW) model |2§] and the Melikhov-Stech (MS) model |25|], both 



based on the constituent quark picture. For the q 2 dependence, the BSW model adopts the 
pole dominance assumption: 

f(i 2 ) = n f ~§ 2W (2-23) 

(1 - q 2 /m 2 ) n 

with m* being the pole mass. The original BSW model assumes a monopole behavior (i.e. 
n = 1) for all the form factors. However, this is not consistent with heavy quark symmetry 
scaling relations for heavy-to-light transitions. The modified BSW model takes the BSW 
model results for the form factors at zero momentum transfer but makes a different ansatz 
for their q 2 dependence, namely, a dipole behavior (i.e. n = 2) is assumed for the form 
factors F 1 , A , A 2 , V, motivated by heavy quark symmetry, and a monopole dependence 
for F Q ,Ai. The experimental value of F® K (0) is about 0.76 |26|], but it is far less certain 



for F^^O). A sum rule analysis |!7[] and in particular a recent lattice calculation JZS] all 
give F® K (0)/F® 7T (0) « 1.20 , in agreement with the results of the BSW and MS models (see 
Table III). 

TABLE III. Form factors in BSW and MS models. 







F^{m 2 K ) 


F 1 DK (m 2 p ) 


F^(m 2 K ,) 


A°K\ml) 




BSW 


0.76 


0.72 


1.01 


1.07 


0.74 


0.77 


MS 


0.78 


0.71 


0.93 


0.91 


0.76 


0.73 



The values of ai and a 2 and their ratio are listed in Table IV. We see that the ratio of 
a 2 /ai is channel dependent, especially its magnitude in TC*tt and Kp decays can be different 
by a factor of 2. However, its phase of order 150° is essentially process independent. To see 
the sensitivity of c^/ai to the IF-exchange contribution, we set £ = and determine 02/^1 
from the isospin analysis of D — » K^n(p) decays. For example, for D — > Ktx decays we 



have 29 



02 

a 1 



V2 



A X/ 2 



13/2 



D^Ktt 



V2 



/2 

h/2 



fn m 2 D - m\ F { 
f K m 2 D -ml Fq 



DK, 



D — *Ktt 



Dtt, 



in 



(2.24) 



K ) 



D^Ktt 
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The results are shown in the last two rows of Table IV. It is interesting to see that for Kit 
and K*tt decays, the phase of 02/01 is about the same as before but the magnitude differs 
slightly, whereas for the Kp system, the magnitude is close to the realistic one but the phase 
is different. This is understandable because the three Kp amplitudes are in phase with one 
another. Hence, 02/01 is real in the absence of the weak annihilation. 

TABLE IV. The parameters a\ and 02 extracted using the BSW model (fist entry) and the MS 
model (second entry) for form factors. Only the central values for the magnitude and the phase 
angle are quoted. To see the sensitivity of 02/01 to the IV-exchange contribution, its value in the 
absence of £ is also shown in last two rows. 





D^Ktt 


D — > K*TT 


D -> Kp 


I«ll 


1.02 


1.23 


0.92 




1.05 


1.28 


0.85 


l°2| 


0.63 


0.53 


0.76 




0.62 


0.45 


0.72 


02/01 


0.62exp(-il52°) 


0.43 exp(-il52°) 


0.82exp(-il48°) 




0.60exp(-il52°) 


0.35exp(-il52°) 


0.85exp(-il48°) 


02/01 (with £ = 0) 


0.88exp(-il49°) 


0.56exp(-il60°) 


-0.87 




0.86exp(-il49°) 


0.46 exp(-il60°) 


-0.90 



What is the relation between the coefficients Oj and the Wilson coefficients in the effective 
Hamiltonian approach ? Under the naive factorization hypothesis, one has 

1 1 

aM = ci(/x) + — c 2 (/i), a 2 (n) = c 2 {p) + — ci(/i), (2.25) 

for decay amplitudes induced by current- current operators Oi,2 (/■*), where Ci^ip) are the 
corresponding Wilson coefficients and N c is the number of colors. In the absence of QCD 
corrections, C\ — 1 and c 2 = 0, and hence class-II modes governed by 02 = l/N c are 
obviously "color-suppressed". However, this naive factorization approach encounters two 
principal difficulties: (i) the coefficients Oj given by Eq. ( |2.25D are renormalization scale and 
75-scheme dependent, and (ii) it fails to describe the color-suppressed class-II decay modes 
due to the smallness of 02- Therefore, it is necessary to take into account nonfactorizable 
corrections: 

Oi = d(ji) + C 2 (/i) + xM^j , «2 = C 2 (/i) + Ci(/i) (J^- + *2(a0) , (2.26) 

where nonfactorizable terms are characterized by the parameters Xi-, which receive corrections 
including vertex corrections, hard spectator interactions involving the spectator quark of the 
heavy meson, and FSI effects from inelastic rescattering, resonance effects, • • •, etc. The 
nonfactorizable terms x«(/-0 wu l compensate the scale and scheme dependence of Wilson 
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coefficients to render a, physical. Using the leading order Wilson coefficients ci(m c ) = 1.274 
and C2(m c ) = —0.529 |3(J for A^g = 215 MeV, where m c (m c ) « 1.3 GeV, it is clear that a 2 
is rather sensitive to \2 and that the nonfactorizable correction to 02 in the Kp system is 
far more important than that in K it decays. 

Empirically, it was found that that the discrepancy between theory and experiment for 
charm decays is greatly improved if Fierz transformed terms in ( |2.25| ) are dropped [p]]. It has 
been argued that this empirical observation is justified in the so-called large- N c approach in 
which a rule of discarding subleading 1/N C terms can be formulated ||. This amounts to 
having universal nonfactorizable terms \i = Xi = —1/N c in Eq. ( [2.26P and hence 



Ql = ci(m c ) « 1.27 , a 2 = c 2 {rh c ) « -0.53 . (2.27) 

This corresponds to a relative phase of 180°. From Table IV we see that the above values 
of ai and a 2 give a good description of the D — > ~K*it decays and differ not too much from 
those values for Kn and Kp systems. Hence, ai and a 2 in the large-iV c approach can be 
considered as the benchmarked values. In the heavy quark limit, nonfactorizable terms Xi 
are calculable due to the suppression of power corrections provided that the emitted meson 
is light while the recoiled meson can be either light or heavy ||. In the QCD factorization 
approach, \i are found to be positive for B decays f| . This means that the empiric large- N c 
approach cannot be generalized to the B system. For charm decays, the charmed quark is 
not heavy enough to apply the QCD factorization approach or pQCD in a reliable manner. 
To our knowledge, the sum-rule approach is more suitable for studying the nonfactorized 
effects in hadronic D decay Q. 

C. Cabibbo-suppressed modes and SU(3) breaking 

So far the diagrammatic amplitudes are determined for Cabibbo-allowed D decays. When 
generalized to Cabibbo-suppressed modes, there exist some sizable SU(3) breaking effects 
which cannot be ignored. Table V shows the predicted branching ratios (see column 3 
denoted by £> t h e0 ryi) for some of the Cabibbo-suppressed modes using the reduced amplitudes 
determined from Cabibbo-allowed decays. By comparing with experiment, we see some large 
discrepancies. To be specific, we consider the ratios: 



Ri 



Vr, 



cd 



T(D+ -> 7T°7r+) T(D° -> K+K-) 



= " fTyn , _ (2-28) 



In the SU(3) limit, R\ — R 2 — 1, while the experimental measurements R\ = 3.39 ± 0.70 
and i?2 = 2.80 ± 0.20 [15] show a large deviation from SU(3) flavor symmetry. 

As first stressed in |JT|, model predictions are very difficult to accommodate the mea- 
sured value of R\. It was originally argued in the same reference that large SU(3) violation 
manifested in R\ can be accounted for by the accumulations of several small SU(3) breaking 
effects, provided that F D7r (0) > F® K (0). However, smaller F D7r (0) is preferred on theoretical 
grounds as discussed before. To accommodate the data of 7r + 7r°, it is clear that one needs 
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TABLE V. Predicted branching ratios of some Cabibbo-suppressed D decays without and with 
SU(3) violation (denoted by the subscripts "theoryl" and "theory2", respectively) and comparison 
with experiment. The amplitude £ q in D° — > K°K° decay denotes the gg-popping VF-exchange 
amplitude. 



Decay mode 


Amplitude 


^theoryl 


^thcory2 


tfcxpt 11 


D+ -> vr+vr 




7.6 x 10" 4 


2.4 X 1(T 3 


(2.5 ±0.7) x 10~ 3 


D° -» 7T+7T- 


y£\T + 6")^ 


2.1 x 10~ 3 


1.7 x i(r 3 


(1.52 ±0.09) x 10~ 3 






1.2 x 10~ 3 


0.9 x 10~ 3 


(8.4 ±2.2) x 10" 4 


D+ -» K+i? U 




4.7 x 1(T 3 


7.4 x 1(T 3 


(7.4 ± 1.0) x 10~ 3 


L>° K+K- 




1.8 x 10~ 3 


4.2 x 10~ 3 


(4.25 ±0.16) x 10~ 3 










(6.5 ± 1.8) x 10~ 4 



T W7r > T and \C W7r \ < \C\. Allowing modest SU(3) violation in the individual quark-graph 
amplitudes, 



% n = 1.25 T, = 1.6 e- il40 ° x 1(T 6 GeV, 

S n7T = 1.9 e il40 ° x 1(T 6 GeV, (2.29) 



the data of D — > titt are well accounted for (see Table V). However, we would like to stress 
that we do not claim that Eq. (2.29) is the solution for Cabibbo-suppressed D — > titt decays; 



we simply wish to illustrate that modest SU(3) violation in each quark-graph amplitudes can 
lead to a large SU(3)-breaking effect for R\. In the factorization approach, however, it is 
difficult to understand why %^ > T. 

The ratio of K + K~ to 7r + 7r~ is a long-standing puzzle. The conventional factorization 
approach leads to R% m 1 when weak annihilation contributions are neglected. In the dia- 
grammatical approach, it is found that R2 can be accommodated by allowing SU(3) violation 
in the tree and ^-exchange amplitudes: 

T KK = 1.25T, £ KK = 1.7 e m ° x 1(T 6 GeV, Akk = A. (2.30) 



It is known that in the limit of SU(3) symmetry, D° — > K°K vanishes. This decay 
receives contributions from inelastic final-state scattering in analog to Fig. 1(a) and it has 
been discussed recently in [|32l,j3~3|l . 



III. WEAK ANNIHILATIONS AND RESONANT FINAL-STATE 

INTERACTIONS 

We learned from Sec. II some important information about the weak annihilation topolo- 
gies S and A: (i) Their contributions are in general comparable to the internal PF-emission 
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topological amplitude and they have a phase of order 90° relative to T with an opposite sign 
between £ and A. (ii) A sign flip of the ^-exchange amplitude may occur in K*n and Kp 
decays and this is very crucial for explaining why D — > Kp amplitudes are in phase with one 
another, but not the case for D — > K*ir and D — » Kir decays. The purpose of this section 
is to explore these two features. 

Under the factorization hypothesis, the factorizable PV-exchange and ^-annihilation am- 
plitudes are proportional to a 2 and a\, respectively. They are suppressed due to the smallness 
of the form factor ^ (m 2 D ) at large q 2 = m 2 D . This corresponds to the so-called helicity 
suppression. At first glance, it appears that the factorizable weak annihilation amplitudes are 
too small to be consistent with experiment at all. However, in the diagrammatic approach 
here, the topological amplitudes C, £, A can receive contributions from the tree amplitude 
T via final-state rescattering, as illustrated in Fig. 1 for D° — > K tt° decay: Fig. 1(a) has 
the same topology as H^-exchange,| while 1(b) mimics the internal py-emission amplitude 
C. Therefore, even if the short- distance PP^-exchange vanishes, a long-distance PP^-exchange 



can be induced via inelastic FSIs |34lj35fl . Historically, it was first pointed out in |J that 
rescattering effects required by unitarity can produce the reaction D° — > K <ft, for example, 
even in the absence of PU-exchange diagram. Then it was shown in J7| that this rescattering 
diagram belongs to the generic PP^-exchange topology. 

There are several different forms of FSIs: elastic scattering and inelastic scattering such 



as quark exchange, resonance formation,- • -, etc. As emphasized in ||34|| , the resonance for- 
mation of FSI via qq resonances is probably the most important one. Indeed, there are two 
indications about the importance of resonant FSIs for weak annihilation topologies: First, 
the sizable magnitude of S and A and their large phases are suggestive of nearby resonance 
effects. Second, an abundant spectrum of resonances is known to exist at energies close to 
the mass of the charmed meson. 

Since FSIs are nonperturbative in nature, in principle it is notoriously difficult to calculate 
their effects. It is customary to evaluate the long-distance PT-exchange contribution, Fig. 



1(a), at the hadron level manifested as Fig. 2 p2| , |36| -p^1 . Take D° — > K°7T° as an illustration. 



Fig. 2(a) shows the resonant amplitude coming from D° — > K~tt + followed by a s-channel 
J p = + particle exchange with the quark content (sd), for example, ^(1950), which 
couples to K°7r° and K~ti + . Fig. 2(b) corresponds to the t-channel contribution with one- 
particle exchange. As discussed before, it is expected that the long-distance PU-exchange 
is dominated by resonant FSIs as shown in Fig. 2(a). However, a direct calculation of 
this diagram is subject to many theoretical uncertainties. For example, the coupling of the 
resonance to Kn states is unknown and the off-shell effects in the chiral loop should be 
properly addressed fl36fl . Nevertheless, as pointed out in |]34] , |39| , most of the properties of 



resonances follow from unitarity alone, without regard to the dynamical mechanism that 



is also pointed out by Close and Lipkin [18] that the prominent weak annihilation may be 
largely due to final-state resonance scattering. 
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(a) 




(b) 

FIG. 1. Contributions to D° — > K°ir from the color- allowed weak decay D° — > ET _ 7r + followed 
by a resonant-like rescattering (a) and quark exchange (b). While (a) has the same topology as the 
VF-exchange graph, (b) mimics the color-suppressed internal VF-emission graph. 

produces the resonance. Consequently, as we shall see below, the effect of resonance-induced 
FSIs [Fig. 2(a)] can be described in a model-independent manner in terms of the mass and 
width of the nearby resonances. 

A. Formalism 

In the presence of resonances, the decay amplitude of the charmed meson D decaying 
into two mesons M\M.2 is modified by rescattering through a multiplet of resonances |P9|p*] 



A(D 



Mi Mj ) resonant- FSI 



— i 



A(D 

r \ ^ (r) (r)* a i 7-j 

E-m R + i T/2^ Y H ( 



M k M h 



(3.1) 



**The same expression for Eq. ( |3.1| ) is also given in [^] except that the phase in Eq. (3.3) of 
is too small by a factor of 2. 
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(a) 



(b) 



FIG. 2. Manifestation of Fig. 1(a) as the long-distance s- and t-channel contributions to the 
VF-exchange amplitude in D° — > i^°7r°. The thick line in (a) represents a resonance. 



where the summation runs over various mass degenerated resonances with the same values 
of the resonance mass uir and width T, and the normalized coupling constants of 

MiMj with the resonance r, obeying the ortho normal relations 



u ij u ij 



M|2 



(3.2) 



The presence of a resonance shows itself in a characteristic behavior of phase shifts near the 
resonance. For each individual resonant state r, there is an eigenstate of A(D 



with eigenvalue [39 



J2i5 r 



l-i 



m D — rriR + iT/2' 



(3.3) 



or 



tan<L 



2(m D -m R ) 



(3.4) 



in the rest frame of the charmed meson. Therefore, resonance-induced coupled-channel 
effects are amenable technically in terms of the physical resonances. 

To illustrate the effect of FSIs in the resonance formation, consider the decays D° — > K^Pj 
as an example. The only nearby + scalar resonance with the (sd) quark content in the charm 
mass region is r = Xq(1950) and the states K^Pj are K~n + ,K tt°,K t],K rf . The quark- 
diagram amplitudes for D° — > K~tt + , K 7r°, K r] q and K 7] s are shown in Table I. It is 
convenient to decompose D — > Kir amplitudes into their isospin amplitudes (see Eq. ( |2.8| ) 
and Table I): 



A(D° -> (Kir) 3/2 ) 



[T + £), A(D° - (Kir) 1/2 ) = —(2T-C + 3£), (3.5) 



where the subscripts 1/2 and 3/2 denote the isospin of the Kit system. Consider the D-type 
coupling for the strong interaction P\P 2 — > P' (P'' scalar meson), namely kTt (P'{P\, P 2 }) 
with k being a flavor-symmetric strong coupling |34|]. Noting that (Kir) 3 / 2 does not couple 
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to (Ktt)i/2, K r) q , and K r/ s via FSIs, the matrix c 2 arising from two D-tjpe couplings in 



the I — \ sector has the form: 



2 2 
C OC K 



in the basis of (Ktt)i/ 2 , K°r] q , ~K°r] s . Since XM c ijl 
matrix c 2 reads 



Y.3 \/3\ 

2 2 V2 \ 

yj i i 

2 2 v 7 ! 

^ -i- 1 / 



/ 3 



(3.6) 



3k 2 , it follows that the normalized 



c 2 ^c 





1 


1 


1 


2 


2^ 


v/6 




1 


1 


1 




2V3 


6 






1 


1 


1 


s/E 


3V2 


3 



(3.7) 



Then it is easily seen that 
i0 —o 



/1(/J U -> X r/ s ) = 5 = e + (e 

x 



2i<5 r 



r 



±=A{D» -> (Ktt) 1/2 ) + ^ A(L>° - A,) + - F\ 



and hence 



£ = e + (e 



2iS r 



t 

6+ 3 



(3i 



(3.9) 



where we have used t, c, e, a to denote the corresponding reduced amplitudes T, C, £, A be- 
fore resonant FSIs. Likewise, it is straightforward to show that T = t and C = c. Therefore, 
resonance-induced FSIs amount to modifying the VF-exchange amplitude and leaving the 
other quark-diagram amplitudes T and C intact. We thus see that even if the short-distance 
VF-exchange vanishes (i.e. e = 0), as commonly asserted, a long-distance ^-exchange con- 
tribution still can be induced from the tree amplitude T via FSIs in resonance formationfl 

Q ' ' 

Likewise, from Cabibbo-allowed decays — > K K + , ir + 7] g and ir + rj s one can show that 
the IF-annihilation amplitude after resonant FSIs reads 



A = a + (e 



2i8 r 



a + 



C 



(3.10) 



^An expression similar to Eq. ( |3.9[ ), 

£ = e + (cos 5e iS -l)(e+ 



was first obtained by Zenczykowski [34| by applying the strong reaction matrix Kq together with the 
unitarity constraint of the S matrix to study the effects of resonance-induced FSIs. However, since 
(e 2lSr — 1) = 2 (cos 5 r e lSr — 1), it is clear that the contribution of resonant FSIs to the VF-exchange 
amplitude as given above is too small by a factor of 2. This has been corrected in |2 
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Note that, contrary to the VT-exchange case, the long-distance VT-annihilation amplitude is 
induced from the color-suppressed internal ly-emission. 

As for the W^-exchange graph in D — > K*n and Kp decays, we consider a 0~ resonance 
P' which couples to VP and PV states with the F-type coupling, k'Tt(P'[V, P]). Proceeding 
as before, the 6x6 normalized coupling matrix reads 



1 







c 


— c 







(3.11) 



in the basis of {K tt)i/ 2 , K T] q , (f>K , (Kp)i/ 2 , K uj and r\ s K , where c is the matrix given 
by Eq. Q and 



A(D° - (K\) 3/2 ) = ^=(T V + S P ), A(D° - (K\) 1/2 ) = -=(2T V - C P + 3£ P 



V3 

A(D° ^(Kp) 3/2 ) = ^=(T P + £ V 



A(D° (Fp) 1/2 ) = ^(2T P - Cy + 3£y). (3.12) 



It is straightforward to show that 7py and Cpy are not affected by resonant FSIs and 



or 



Si 



e P + -le 



2i5 r 



-1) 



Sv = e v + ^(e 2i5r - I) 



e P -e v + -(T v - T P ) 
e v -e P + ^(T P - T v ) 



(3.13) 



E P + E v = e P + ey, 



£p-£v = ep-e v + {e 2tdr -l)(e P -e v - -(T P - T v ] 



(3.14) 



Note that the terms in square brackets in Eq. (|3.13| ) have opposite signs for £p and Ey owing 
to the charge-conjugation of the strong coupling.^ We will come back to this point later. 

As for the ^-annihilation amplitudes Ap and Ay, a direct analysis of resonant FSIs in 
Cabibbo-allowed decays D+ — > pn, p + rj(r]'), cj7r + ,07r + , K K shows that the reduced ampli- 
tudes Tpy and Ap + Ay are not affected by FSIs in resonance formation 



Tp — t P , 



Ty — ty, 



Ap + Ay = ap + ay. 



(3.15) 



The relevant normalized coupling matrix in the basis of (pvr)i, K K + and K* + K is given 
by 



/ 2 I — I \ 

3 3 3 1 

1 1 _ 1 

3 6 6 

_ 1 _1 1 

3 6 6 



(3.16) 



"The expression of Eq. ( |3.13| ) or ( |3.14| ) differs from the results obtained in p4|J?2f for £p and By. 
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where 



A(Dj -> Wi) = - A + ) - A(£>+ - P + A\ =A P - A v , 



- ( P tt) 2 ) = ^ A + ) + - P + «°)] = 0- (3-17) 



It follows that 

Ap + Ay = ap + a v , 



A r - A v or - «r {e 2i5r -l)(a P -a v + -{C P -C v ) ). (3.18) 



and Cpy are not affected. This result was first obtained by Zenczykowski ||34j| . Note the 
similarity between the expressions of (|3.18|) and (|3.14|). 



B. Phenomenological implications 

1. Resonance-induced weak annihilations 



Eqs. (|3.9|) , ( p. 10|) , (|3.14|) and (|3.18|) are the main results for weak annihilation amplitudes 



induced from FSIs in resonance formation. We see that even in the absence of the short- 
distance weak annihilation, a long-distance weak annihilation can be induced via resonant 
FSIs. For parity- violating D — > PP decays, there is a J p = + resonance i^o(1950) in the sd 
quark content with mass 1945 ± 10 ± 20 MeV and width 201 ± 34 ± 79 MeV [pj. Assuming 
e = in Eq. ( |3.9D , we obtain 

8 = 1.43 x 10~ 6 exp(zl43°) GeV, (3.19) 

which is close to the "experimental" value S = (1.60 ±0.13) x 10 -6 exp(ill4°) GeV [cf. Eq. 
( |2.9|) ]. Presumably, a non- vanishing short-distance e will bring the phase of £ in agreement 
with experiment. Resonance-induced coupled-channel effects are governed by the width and 
mass of nearby resonances, which are unfortunately not well determined. For example, a 
reanalysis in a .fT-matrix formalism |^0[ quotes nip = 1820 ± 40 MeV and T = 250 ± 100 



MeV for the same resonance. This leads to £ = 1.67 x 10 6 exp(— zl58°) GeV. Therefore, 
we conclude that one needs a + resonance heavier than the charmed meson. Contributions 



from more distance resonance at 1430 MeV are smaller (see also |38|.|TT[ ) . 

Since a nearby + resonance cto in the charm mass region has not been observed, we 
will not make an estimate of the resonance-induced ^^-annihilation amplitude in Cabibbo- 
allowed D — > PP decay. In the factorization approach, it is expected that £/A ~ 02/^1 
and hence |^4| > \£\. Nevertheless, it is clear from Eq. (|3.10 ) that the long-distance W- 



annihilation is slightly smaller than the W-exchange one because the former is induced from 
the color-suppressed amplitude C and this is consistent with Eq. 
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To estimate the resonance- induced W- annihilation amplitude Ap — Ay, we employ 
7r(1800) as the appropriate 0~ resonance with rap = 1801 ± 13 MeV and T = 210 ± 15 
MeV JL5|. Assuming ap — ay = 0, we find from Eq. (|3.18| ) that 

\A P -A v \ = 2.4 x 10^ 7 , (3.20) 



which is in agreement with the experimental value given by Eq. ( 2.18 ). Note that if we set 
ap = ay = 0, this will lead to Ay = —Ap which will imply a vanishing — > tt + uj, in 
contradiction to the experimental observation. 



2. Hadronic D decays to rj or if 



Weak annihilation effects are crucial for some two-body decays involving one single isospin 
component, e.g. the final state containing an 77 and vj '. This has been discussed in detail 
in |2^j2^j4T|]. To see this, we consider the Cabibbo-allowed decays D° — ► "K°(r],'r)') and 
Cabibbo-suppressed modes D + — > ir + (j],r)'). In realistic calculations we use the 77—7/ mixing 
angle 9 = —15.4°, but it suffices for our purposes to use 9 = —19.5° to discuss the essence of 
physics. The quark-graph amplitudes then read 



A(D C 



A(D H 



Wrf) = C/V3, 
2C + 2A), 



A(D° - 
A(D + 



K°r]>) 

■ 7T + V) 



(C + 3£)/V6, 

Yl^(X-C + 2A). (3.21) 



Since £ is comparable to the color-suppressed C in magnitude, the decay D° — > IT'rj' is 
largely enhanced by jy-exchange. We see from Table VI that its branching ratio is enhanced 
by resonance- induced FSIs by almost one order of magnitude, whereas D° — > K°r] remains 
essentially unaffected. Therefore, we conclude that it is the ^-exchange effect that accounts 
for the bulk of B(D° — > K 77') and explains why K 7/ > K rj. 

As for D + — > 7t + (t7, rj') decays, it is clear from Eq. (|2.9| ) that the interference is construc- 
tive between 2A and T + 2C and destructive between 2A and T — C. Hence, the presence of 
ly-annihilation is crucial for understanding the data of D + — ► tt + 7] (see Table VI). 

TABLE VI. Branching ratios (in units of %) of the charmed meson decays to an 77 or rj ' . 



Decay 

w 

D° 
D+ 



K rj 
K rj 

> 7T + ?7 



D+ -> ir+rj' 



Without weak annihilation 
0.64 
0.31 
0.09 
0.27 



With weak annihilation Expt. [15] 

0.66 0.70 ±0.10 

1.76 1.71 ±0.26 

0.37 0.30 ± 0.06 

0.39 0.50 ±0.10 
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C. Difficulties with weak annihilation amplitudes Spy and Apy 



For the VT-exchange amplitude in D — > VP decays, at first sight it appears that its sign 
flip from £p to Ey can be naturally explained since the 0~ resonance is expected to have equal 
and opposite couplings to K*~n + and K~ p + , as shown in Eq. ( |3.13|) . It is obvious from Eq. 
(|3.14|) that Sp = —Sy in the absence of short-distance VF-exchange contributions. However, 



a further study shows some problems. First, a possible candidate of the resonance near 
m D is the #(1830) with mass ~ 1830 MeV and width ~ 250 MeV [|15j, but only the K<\> 
decay mode has been reported. Second, due to the F-type coupling of the resonance with 
VP states, the resonance contribution from the leading tree amplitude T is proportional to 
the difference Tp — Ty, which is small. Assuming ep = ey = for the moment, we obtain 
from Eqs. (|2~T0| ), (|2~lp and d3~l3|) that 



S P = -S v = 1.6 x 10~ 7 exp(il7°) (3.22) 

from the resonance #(1830). Evidently, Sp and Sy are not governed by resonant FSIs, 
contrary to the original conjecture advocated in |1(J. In order to explain the sign flip of the 
H^-exchange amplitude, it appears that the (short-distance) V^-exchange amplitudes ep and 
ey have to be sizable in magnitude and opposite in signs or the 0~ resonance couples strongly 
to one of the K*~~tt + and K~ p + states, or the long-distance t-channel effect analogous to 
Fig. 2(b), which has been ignored thus far, gives the dominant contributions to Sp and Sy 
with opposite signs. In any of these cases, it is not clear what is the underlying physics. 
Therefore, the sign flip of the ^-exchange amplitude in Cabibbo-allowed VP decays remains 
unexplained. 

We also face some difficulties for understanding the iy-annihilation amplitudes Ap and 
Ay in D — > VP decays. This is because naively one will expect — ► n + u) is suppressed rel- 
ative to .D+ — > 7i + p°. The argument goes as follows. The direct H^-annihilation contributions 
via cs —>■ W —>■ ud are not allowed in — > tt + u), p + r], p + rj' decays since the (ud) has zero 
total angular momentum and hence it has the quantum number of tt + . Therefore, G(ud) = — 
and the final states should have odd G-parity. Since G-parity is even for U7i + and odd for 
7T + p°, it is clear that the former does not receive direct H^-exchange contribution. Can one 
induce — > h + oj from resonant FSIs ? The answer is no because the J = 0, 1 = 1 meson 
resonance made from a quark-antiquark pair ud has odd G parity. As stressed in |42||, the 



even-G state tt + uj (also prj and prj') does not couple to any single meson resonances, nor to 
the state produced by the VT-annihilation diagram with no gluons emitted by the initial state 
before annihilation. This is indeed consistent with Eq. (|3.18|) which states that, contrary to 
Ap — Ay, Ap + Ay does not receive any qq' resonance contributions. Therefore, it will be 
expected that D+ — > n + u is prohibited (or Ay ~ — Ap), whereas — > 7r + p° receives both 
factorizable and resonance-induced iy-annihilation contributions. Experimentally, it is the 



other way around: B(D+ -> n + u) = (2.8 ± 1.1) x 10" 3 |4| and B{D+ -> n+p°) ~ 6 x 10" 4 



19| . Hence, it is important to understand the origin of the M^-annihilation contribution. As 



noted in passing, Eqs. ( |2.17|) and ( |2.18|) suggest that the phase difference between Ap and 
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Ay is less than 90° and that the magnitude of Ap or Ay is smaller than Sp and Sy, contrary 
to the PP case. The suppression of PF-annihilation is expected because of the G-parity 
constraint. Since ^-annihilation occurs only in the system for Cabibbo-allowed decays, 
its suppression may help to explain why t(D£) > t(D°) at the two-body or quasi-two-body 
decay level. 



IV. COLOR-SUPPRESSED INTERNAL ^-EMISSION AMPLITUDE 



The relative phase between T and C indicates some final-state interactions responsible for 
this. Fig. 1(b) shows that final-state rescattering via quark exchange has the same topology 
as the color-suppressed internal VF-emission amplitude. At the hadron level, Fig. 1(b) is 
manifested as FSIs with one particle exchange in the t channel |]44]j32| , see Fig. 3. Note that 
although Fig. 3 is very similar to Fig. 2(b), the exchanged particles here are K* and K 
rather than p and ir. Another approach is based on the Regge pole model [02| . Admittedly, 
the estimate of Fig. 3 is subject to some theoretical uncertainties as discussed before, e.g. 
the off-shell form factor appearing in the chiral loop calculation is unknown. Nevertheless, 
we still can learn something useful. For example, Li and Zou JH] have calculated rescattering 
FSIs for D+ -> K ir+ and D+ -> K p + via single pion exchange. The absorptive part of 
the chiral loop diagram in analog to Fig. 3 gives the imaginary contribution to the FSI 
amplitude. It was found the imaginary part is negative for Cp and Cy (see Table II of P4]). 
This lends a support to the solution set (i) in Eq. ( |2.16| ). We will leave a detailed study of 
this inelastic FSI effect to a separate publication. 




(b) 



FIG. 3. Manifestation of Fig. 1(b) as the long-distance i-channel contributions to the 
color-suppressed internal VF-emission amplitude in D° — ► K°tt°. 



It is often argued that the relative phase between T and C can be understood in the 
factorization approach as arising from the rescattering through modifications of the phases 
of isospin amplitudes Hl4 , ^5 ,10|. First, one identifies the isospin amplitudes with the fac- 



torizable amplitudes before the final-state phases are turned on in Eq. ( |2.8|) . For example, 
A{D + -> K 7T+) = ^^3/2 = T f + C f . This leads to the isospin amplitudes 

A1/2 = ^(2Z> -C f + 3E f ), A 3/2 = -*=(2) + C f ), (4.1) 
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where Tj,Cf,Ej are the factorizable tree, color-suppressed, and VT-exchange amplitudes, 
respectively. Second, after introducing isospin phases to Eq. ( J2.8| ), it is straightforward to 
show that 



T + S = (T f + E f )e l5 ^ --(T f + C f )(e iS ^ 2 - e^ 2 ), 

C-S = (C f - E f )e iSl / 2 ~-{T f + C f )(e iS V 2 - e^ 3 / 2 ). 

3 



(4.2) 



Therefore, the different phases of C and T are a consequence of rescattering. This has an 
important implication for the class-II mode D° — > K°-k°. Even if the factorizable Cf and 
Ef amplitudes are small, the weak decay D° — > K~n + followed by the inelastic rescattering 
K~n + — > ~K°7t° can raise B(D° — > ~K°it ) dramatically and lower B(D° — > K~tt + ) slightly. 

Of course, the above picture for rescattering is too simplified and it does not offer a 
genuine explanation of the phases of quark-graph amplitudes C and T . First, the decays 
~K°7] and K°r]' have only one isospin component and the above analysis does not give a clue 
as to why they are not color suppressed. Second, the isospin phase difference in D — > Kp 
decays is near zero and yet a relative phase of order 150° between C and T is found in the 
diagrammatic approach [see Eq. ( 2.16Q 1- Third, the above isospin analysis has no power of 
prediction as the phase difference is unknown from the outset; even if elastic Kir scattering 
is measured at energies y/s = m D , the isospin phases appearing in ( |2.8|) and ( |4.2|) cannot be 
identified with the measured strong phases.0 



V. COMPARISON WITH B DECAYS 



It is instructive to compare the present study with the B decays. To proceed, we quote 
some of the results for B — > Dir, D*tt decays f29[ 



B -»• D 



7T 



B -> D*ti 



Ai/2 

A3/2 

Al/2 



A 



.3/2 



1.00 ±0.14) e i29 °, 
1.05 ±0.10) e i29 °, 



^ = (0.40~0.65)e i59 °, 



a 2 



(0.30 ~ 0.55)e 



«63° 



(5.1) 



The relative phase between a\ and a 2 is of order 60°. The QCD factorization approach 
§ implies that 5 1/2 - 5 3/2 _ = 0(A QC D/m Q ) and l 1/2 /(v / 21 3 / 2 ) = 1 + 0(A QC D/m Q ). We 
see that, except for D — > Kp, the isospin phase difference indeed decreases from charm (of 
order 150°) to the bottom system and the ratio of isospin amplitudes in D decays shows 



§lf there are only a few channels are open as the case of two-body nonleptonic decays of kaons and 
hyperons, the isospin phases there (or decay amplitude phases) are related to strong-interaction 
eigenphases. However, when there are many channels open and some channels coupled, as in D 
and especially B decays, the decay phase is no longer the same as the eigenphase in the S'-matrix. 
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a sizable departure from the heavy quark limit. The relative phase of a% and a 2 is crucial 
for understanding the destructive interference in the class-Ill mode D + — > K°tt + and the 
constructive interference in B~ — > D°7r~ . 

Since the color-suppressed mode B — > J/ipK does not receive any weak annihilation 
contribution and the penguin contribution to this decay is very tiny, the parameter a 2 can 



be directly determined from experiment with the result \a 2 (B — > J/ipK)\ = 0.26 ±0.02 [|5 
It is evident that even in B decays, a 2 varies from channel to channel. 

In obtaining Eq. ( |5.1|) we have neglected the H^-exchange amplitudes in B — > D^n 
decays. It is generally argued that weak annihilation is negligible as helicity suppression 
should be more effective due to the large energy release in B decays. Owing to the absence 
of nearby resonances in the B mass region, the weak annihilation amplitude will not be 
dominated by resonance-induced FSIs, contrary to the charm case. Hence, the formalism 
developed in Sec. Ill for D decays is not applicable to B mesons. 

Recently, there are some growing hints that penguin-induced weak annihilation is impor- 
tant. For example, the theory predictions of the charmless B decays to K*tt, Kp and K* r\ 
based on QCD factorization are too small by one order of magnitude in decay rates |[46|| . 
This implies that the weak annihilation (denoted by P e or P a in the literature) may play 
an essential role. Indeed, a recent calculation based on the hard scattering pQCD approach 
indicates that the K*tt modes are dominated by penguin-induced annihilation diagrams |[47|| . 
This is understandable because the usual helicity suppression argument works only for the 
weak annihilation diagram produced from (V — A)(V — A) operators. However, weak anni- 
hilations induced by the (S — P)(S + P) penguin operators are no longer subject to helicity 
suppression and hence can be sizable. 

It is anticipated that the soft FSI contributions to the color-suppressed topology C are 
dominated by inelastic rescattering [48]. As for the phase of the ratio of a 2 /ai, the rescatter- 
ing contribution via quark exchange, D + tt~ — > D°ti , to the topology C in 5° — > D°ti has 
been estimated in using p trajectory Regge exchange. It was found that the additional 
contribution to D°rr from rescattering is mainly imaginary: a 2 (D11) / a 2 (Dn) without fsis = 
1 + 0.61exp(73°). This analysis suggests that the rescattering amplitude can bring a large 
phase to a 2 {Dn) as expected. 



VI. CONCLUSIONS 

We have presented a study of hadronic charm decays within the framework of the dia- 
grammatic approach. We draw some conclusions from the analysis: 

1. Based on SU(3) symmetry, many of the topological quark-graph amplitudes for 
Cabibbo-allowed D decays can be extracted from the data. The ratio of a 2 ja\ has 
the magnitude of order 0.60, 0.40 and 0.83 for D — > Kir, K*ir, Kp decays, respec- 
tively, with a phase of order 150°. This implies that nonfactorizable corrections to Kp 
are far more important than K n. 
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2. Except for the H^-annihilation topology in VP decays, the weak annihilation (W- 
exchange or W^-annihilation) amplitude has a sizable magnitude comparable to the 
color-suppressed internal I^-emission with a large phase relative to the tree ampli- 
tude. It receives long-distance contributions from nearby resonance via inelastic final- 
state interactions from the leading tree or color-suppressed amplitude. The effects 
of resonance-induced FSIs can be described in a model independent manner and are 
governed by the mass and decay width of the nearby resonances. 

3. Weak annihilation topologies in D — > PP decays are dominated by nearby scalar 
resonances via final-state resacttering. In contrast, W-exchange in VP systems receives 
little contributions from resonant final-state interactions. 

4. The experimental data indicate that the three decay amplitudes of D — > Kp are essen- 
tial in phase with one another. This requires that either the W-exchange amplitudes 
in Kp and K*ir have opposite signs or the relative phase between the tree and color- 
suppressed amplitudes flips the sign. While the latter possibility is probably ruled 
out by the measurement of D° — > ~K*°ri and by the model calculation of the phase 
of the color-suppressed amplitude, the first possibility is hampered by the observation 
that dominance of nearby resonances is not operative for the W-exchange contribution 
in Cabibbo-allowed D — > VP decays. Therefore, why the sign of the ^-exchange 
amplitude flips in D — > K*ir and Kp decays remains unexplained. 

5. Owing to the G-parity constraint, the W- annihilation amplitude Ap or Ay in D — > VP 
decays is suppressed relative to W^-exchange, contrary to the D — > PP case where W- 
exchange and W^-annihilation are comparable. Since PF-annihilation occurs only in the 

system for Cabibbo-allowed decays, this may help to explain the longer lifetime of 
D+ than D°. 

6. Weak annihilation contributions are crucial for understanding the data of D° — > K°i]' 
and D + — > n + r]. 

7. The relative phase between the tree and color-suppressed amplitudes arises from the 
final-state rescattering via quark exchange. This can be evaluated by considering the 
t-channel chiral-loop effect or by applying the Regge pole method. Much more work 
along this line is needed. 

8. Some Cabibbo-suppressed modes exhibit huge SU(3)-flavor symmetry breaking effects. 
It can be accounted for by the accumulation of several modest SU(3) violations in 
individual quark-graph amplitudes. 
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